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Abstract 

Thouless’s quantum adiabatic pumping is of fundamental interest to condensed-matter physics. 
It originally considered a zero-temperature equilibrium state uniformly occupying all the bands 
below a Fermi surface. In the light of recent direct simulations of Thouless’s concept in cold-atom 
systems, this work investigates the dynamics of quantum adiabatic pumping subject to dephasing, 
for rather general initial states with nonuniform populations and possibly interband coherence. Us¬ 
ing a theory based on pure-dephasing Lindblad evolution, we find that the pumping is contributed 
by two parts of different nature, a dephasing-modified geometric part weighted by initial Bloch 
state populations, and an interband-coherence-induced part compromised by dephasing, both of 
them being independent of the pumping time scale. The overall pumping reflects an interplay of 
the band topology, initial state populations, initial state coherence, and dephasing. Theoretical 
results are carefully checked in a Chern insulator model coupled to a pure-dephasing environment, 
providing a useful starting point to understand and coherently control quantum adiabatic pumping 
in general situations. 

PACS numbers: 03.65.Yz, 73.43.-f, 05.30.Rt, 32.80.Qk 
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Introduction. Thouless’s seminal work on adiabatic quantum pumping meg establishes 
a deep connection between band topology and quantum transport. In particular, the sum 
of the Chern numbers of all the filled bands below the Fermi surface determines the number 
of pumped charges over one adiabatic cycle in a one-dimensional (ID) periodic lattice. To 
date various quantum pumps have been proposed to study topological phases and topo¬ 
logical phase transitions [MS]. Thoulcss’s concept, which can be deemed as a dynamical 
version of the integer quantum Hall effect, is directly simulated this year in two cold-atom 
experiments [131 H]- Thouless’s pumping has also been extended to periodically driven 
quantum systems to manifest the topology of Floquet quasi-energy bands instead of energy 
bands HMT]. 

In simulating quantum pumping with cold-atom mm and waveguide [4] systems, the 
initial state might not be the zero-temperature equilibrium state considered by Thouless. 
There are at least three reasons why nonequilibrium initial states should be investigated 
theoretically. First, if a bosonic system is pumped, there does not exist a Fermi surface to 
automatically guarantee full and uniform band filling. Second, even for a fermionic system, 
loading the particles into a lattice in an actual experiment is nontrivial and might not 
result in uniform populations within one band H3j. Third, a simple initial state localized 
in an optical or waveguide lattice may possess interband coherence (IBC) [1TJ in the band 
representation. These nonequilibrium situations occur even more frequently in studies of 
quantum adiabatic pumping in periodically driven systems It is hence of both 

fundamental and practical interest to extend Thouless’s pumping to nonequilibrium initial 
states. Indeed, in a recent study, IBC is found to induce a remarkable correction to adiabatic 
pumping in periodically driven systems D2|. 

In this Letter we study dephasing effects on adiabatic pumping dynamics with rather 
general initial states. On the one hand, realistic systems are always subject to dephasing, 
so how dephasing affects quantum pumping is of theoretical interest. On the other hand, in 
well controlled experimental systems, it is possible to deliberately introduce dephasing effects 
into the setup so as to manifest the findings described below. Specifically, with a minimal 
pure-dephasing model and certain assumptions, we find that the number of pumped parti¬ 
cles over one cycle comprises two components of different nature: one dephasing-modified 
geometric component weighted by initial populations on different Bloch states, and a second 
component induced by IBC but compromised by dephasing, with both of them indepen- 


2 


dent of the pumping time scale. The overall pumping hence depends on the band topology, 
initial Bloch state populations, initial IBC and dephasing, offering a stimulating starting 
point to understand quantum adiabatic pumping with dephasing and in nonequilibrium sit¬ 
uations. Our theory is fully checked in a Chern insulator model coupled to a pure-dephasing 
environment. 

Dephasing-induced correction to nonadiabatic population transfer. Our theory starts from 
a slowly varying Hamiltonian H(s), where s = vt (it is straightforward to extend to other 
time dependence of s) represents a tunable system parameter, t is the time variable, and v is 
the sweeping rate of s. All variables are assumed to be scaled and dimensionless. The Planck 
constant h is set to unity throughout. For convenience s = vt is also used to reflect the time 
t = s/v when appropriate. An adiabatic protocol starts at s = so and ends at s = s\. We 
assume that for any value of s 6 [so, si], the spectrum of H(s) is never degenerate. Hence 
the conventional quantum adiabatic theorem applies when the sweeping rate v —> 0. For 
a small but finite value of v, the lowest-order nonadiabatic corrections should be carefully 
studied in order to capture the main physics of adiabatic pumping. 

To account for dephasing in a solvable manner in the context of adiabatic pump¬ 
ing, we exploit the following master equation in the Lindblad form [LB]: v4-p(s) = 

— i[H(s)p(s) — p(s)H(s )] + 7 |A(s)p(s)A(s) — \ [A 2 (s)p(s) + p(s)A 2 (s)]j. Here s is already 
used to reflect the time, p(s) is the density matrix of the system at t = s/v, A(s) is a Her- 
mitian Lindblad operator, and 7 denotes the dephasing rate. Though this treatment does 
not explicitly consider bath degrees of freedom and has ignored all non-Markovian effects, it 
does allow us to closely follow our previous work HU and to clearly reveal the competition 
between IBC and dephasing. Further, to have a simplest model incorporating dephasing 
effects, the time-varying Hamiltonian H(s) and the Lindblad operator A(s) are assumed 
to be always commutable. By this construction, the system-environment coupling does not 
directly introduce transitions between instantaneous eigenstates of H(s). That is, we adopt 
a model based on pure-dephasing Lindblad evolution [19]. The instantaneous eigenstates of 
H(s) are denoted by |m(s)), with H(s)\m(s)) = E m (s)\m(s)). Because A(s) commutes with 
H(s), |m(s)) are also chosen as eigenstates of A(s), with A(s)|m(s)) = A m (s)\m(s)). 

We first aim to show how both nonadiabaticity and dephasing jointly change the pop¬ 
ulations on, and the coherence between, instantaneous eigenstates of H(s), up to the first 
order in v during a pumping protocol. We project p(s) onto the instantaneous eigenstates 
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(m(s)| and \j(s)) of H(s), yielding p(s) = Pmj(s)\m(s))(j(s)\. This time-dependent 

representation for system’s density matrix will be used throughout. To the first order in v, 
we obtain the off-diagonal elements (coherence) (see Appendix [2D]) 


— s± m ys) - ’ 

as well as the diagonal elements (populations) 



( 1 ) 


Pjj( s i) Pjj( s o) 



+ c.c. 


-2v E \.Pjj( s o) Pmm («o)] / B (s) ds 


( 2 ) 


m/j 


where B(s) 



■ f]mj (s) — T/ m (s) T/j(s), qw/^) 


| [A m (s) — Aj(s)] 2 . ft is enlightening to discuss the nonadiabatic transition probabilities in 
Eq. There it is seen that A Pj = Pjfis i) — Pjj(s 0 ) is composed of two parts: one part 


related to initial state coherence p m j( s o), an d a second part depending on the differences in 


the initial populations on eigenstates of H(sq). The coherence properties of the final density 
matrix do not appear because p m j(si)(m ^ j), the off-diagonal terms of p(s i), have been 


washed off by dephasing. Interestingly, the initial-state-coherence induced correction to A Pj 
does persist in the presence of dephasing, with the gap function g m j(0) however replaced by 
g m j( 0) — *7mj(0). Furthermore, the second part of A Pj shows that, even without initial state 
coherence, dephasing alone can also induce a nonadiabatic correction proportional to v. As 
such, both initial state coherence and dephasing can affect the proportionality coefficient 
A Pj/v, a clear interplay expected to influence strongly on adiabatic pumping. Note also 
that if only one single eigenstate of H(so) is populated at s = so, then the above result 
is fully consistent with the one derived by Avron et al for a pure-dephasing Landau-Zener 
model [ID] . Numerical results in two- and three-level systems fully confirm the theoretical 
results here m- 

Dephasing-modified adiabatic pumping. In adiabatic pumping, nonadiabatic corrections 
of the order v could accumulate over a pumping cycle, potentially yielding an overall outcome 
independent of v or the pumping time scale m- Consider then a particle moving in a ID 
periodic lattice, subject to pure-dephasing Lindblad evolution. The Hamiltonian is assumed 
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to be H k (s), with k 6 [— tt, tt) being the quasimomentum and H k (s ) = H k (s + 2tt). The 
parameter space formed by k and s is hence a 2-dimensional (2D) torus. The eigenvalues 
of H k (s) are assumed to form well-gapped Bloch bands. A pumping cycle can be realized 
by sweeping s slowly from so = 0 to Si = 2tt with a speed v. The Lindblad dephasing 
operators for the fc-component of the system are assumed to be A k (s), which preserves the 
translational invariance. Due to this treatment, k is conserved during the pumping and 
states with different k are always independent of each other. Therefore, Eqs. (|TJ) and ([2]) 
are applicable for each individual k value. As elaborated in Appendix, the spectrum of 
H k (s ) and A k (s), as well as the initial state populations of the system, are all assumed to 
have reflection symmetries in the fc-space so as to highlight all the pumping terms that are 
independent of v. 

The number of particles pumped through per unit cross section at the end of the pumping 
cycle can be written as a 2D integral, namely, Q = A dk ds f ks , with 


fka = Tr [u 1 p(s)v k (s)] , (3) 

where v k (s) = dH f k is the group velocity operator of a particle with quasimomentum k. 
To evaluate f ks , both the diagonal and off-diagonal elements Pjj(s) and p m j(s) are needed. 
Indeed, f ks = v~ l E, Pjj( s ){j( s )\ v k(s)\j(s)) +u _1 P™j( s )(j( s )\ v k(s)\m(s)) . A straight¬ 

forward but somewhat tedious application of Eqs. (jl| and @ ra then yields (to the Erst 
order of v) four subterms of f ks , i.e., f ks = fj£ + fj® + fj£ + fj? s \ with 


- mj (s)R.e 

c(a) _ - ' 


fl7 = 2j>,(0)- 


t dj(s) 
^ dk 


m{s))(m(s)\^) 


m+j 


rvw +1 


/l? = 25 >,(<v 


Im 




dfc 


3,m 

m^j 


r 2 ml (s) +1 


/£’ = -^2^dRe[p mi (0)C( S = 0)] 


j,m 

m^j 


Ad) _ tnl ^'( 0 )| dj dfc 0) | m ( 0 ^ 

hs 7mj (0) + ig m j (0) ’ 


J,m 

m.^j 


( 4 ) 


where C(s) = 




9mj(s)[7rnj(s)+igmj(s)]l 


's) = 7mj , and all other quantities do carry a k- 

" 7 9mj \ s ) 
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dependence but not spelled out explicitly, For example, p m j(s) now means the density matrix 
element already projected onto the subspace with quasimomentum k, \j(s)) above refers to 
an eigenstate of H k (s), T mj (s) should also be understood as a function of k, with 7 m j(s) 
and g m j(s) defined at the same individual k values. As seen above, all the four subterms 
of fks are now independent of v. That is, accumulating the transport behavior over one 
entire pumping cycle allows us to capture all the subtle nonadiabatic and dephasing effects 
proportional to v. Equation Q goes beyond a previous result [2TJ because we not only 
incorporate interband dephasing (as reflected by T m j) as a function of the band gaps, but 
also account for IBC. 

Each of the four subterms in Eq. Q should be discussed. We do so mainly by asking 
what happens to them if removing dephasing [7mi(s) = 0, T m j(s) = 0]. First, depending 
on initial populations would vanish, hence a subterm entirely due to dephasing. Second, 
expectedly reduces to an integral over Berry curvatures of Bloch wavefunctions weighted 
by their initial populations. Third, fjfj would recover an expression parallel to that in 
Ref. [FT], where how IBC in driven systems corrects adiabatic pumping was first studied. 
This is especially encouraging because our derivation here refers to a nondriven system and 
also takes an entirely different route than Ref. HZ]. The alternative derivation in Ref. na 
clearly indicates that f ks is of a dynamical, not a geometrical origin. The found IBC effect is 
seen to persist well in the presence of dephasing, insofar as its explicit form is only modified 
by dephasing via g m j(s) —* g m j(s)—i^ m j(s) [see the expression of C(s) defined below Eq. Q]. 
Lastly, in obtaining flf we already assumed that the off-diagonal elements p m j{si = 2 tt) 
have decayed to zero m- Hence naively setting 7 mj (0) = 0 alone in the expression for 
flf does not suffice to recover any unitary limit. Because this last subterm does not even 
depend on the pumping protocol (no derivatives with respect to s is involved), this subterm 
largely originates from a current inherent in the initial state itself, whose contribution to 
the pumping is accumulated only for a dephasing time scale. Overall, the four subterms 
presented in Eq. Q can now be used to predict, both qualitatively and quantitatively, the 
features of adiabatic transport emanating from rather general nonequilibrium states. 

Adiabatic pumping in a Chern-insulator model. We use the Qi-Wu-Zhang model Hamil¬ 
tonian [ 22 ] . 

H k (s ) = sin(fc)cr x + sin(s)cTj / + [5 + cos (k) + cos(s)]cr 2 , (5) 

which contains the main feature of Chern insulators [23j. This system can describe spin- 
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1/2 ( cr x ,y jZ are Pauli matrices) fermions with spin-dependent nearest-neighbor hoppings on a 
square lattice, with 6 being an energy bias parameter. Though originally k and s refer to two 
quasimomenta along different directions, here we may also view the above Hamiltonian as a 
ID system with s being an experimentally tunable system parameter. The instantanenous 
eigenstates of Hk(s) are represented as an d |'0^(s)), with eigenvalues E k ,i(s) < 

Ek, 2 (s). To introduce pure dephasing, we assume the Lindbald operators A k (s) to be the 
same as H k (s). The spectrum of H k {s) [hence also of Afc(s)] is indeed symmetric in k. 
Remarkably, Eq. Q applied to this two-band model yields particularly simple expressions 
for fj£ and f£\ with 


fks = A Pk( 0 ) 


21 


r 2! + l 


G 


ks j 


fS = Am(0)=5— -Sl k „ 


q, + i 


( 6 ) 


where T 2 i, as defined before, is now given by T 2 i = 7 [Efc i 2 (s) — Ej^s)] /2. Apfc(0) in Eq. ([ 6 ]) 
refers to the initial population difference between the ground and excited bands. Here 


Gks = 2Re 


and n ks = 21 m 




are 


[{ l k v4 1} ( s ) iv4 2) (*)> (?/4 2) (s) 1 ^v4 1} ( s )) 

the Fubini-Study metric of the ground state bundle and the more familiar Bloch band Berry 
curvature [21J. Clearly then, without IBC, the adiabatic pumping solely determined by 
and ffe is of geometrical nature, but modified by dephasing. The whole story of adiabatic 
transport under a pumping protocol is however completed by fjfj and flf, whose explicit 
expressions follow Eq. (|4]) and are not given here. 

To verify our theory, we numerically evolve various initial states at each quasimomentum 
value using the above-mentioned pure-dephasing Lindblad master equation. We then obtain 
Q by integrating the numerically found f ks over k and s. Consider Erst an initial state 
uniformly filling the bottom band only. I 11 this case there is no IBC, so the pumping is 
entirely determined by and ff® in Eq. ( 6 ). As shown in Fig. [lj theory and numerics 
are in perfect agreement. For the topologically nontrivial case with S = 1 (top panel), Q 
decreases from a quantized Chern number (Q = 1) to zero as the dephasing rate 7 increases. 
Thus, the topological relevance to adiabatic pumping is gradually suppressed by dephasing. 
For the topological trivial case with S = 2.5 (bottom panel), however, Q first decreases 
and then increases with 7. We stress that this dephasing-induced pumping is still directly 
connected with the G ks metric. The non-monotonous dependence of Q on the dephasing 
rate 7 is easily understandable from the expression of ft ^ 1 , which reflects a competition 
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FIG. 1: (color online) Number of pumped particles Q vs. the dephasing rate 7 for (a). 6 = 1.0, 
(b). <5 = 2.5. The system Hamiltonian is described in Eq. ([ 5 ]). The adiabatic sweeping rate is 
assumed to be v = 1 CF 3 in our numerical calculations. 


between the dephasing-modified (effective) band gaps and dephasing-modified nonadiabatic 
transition rates. 

We next consider initial states that coherently populate the two bands. As the first exam¬ 
ple, the initial density matrix is chosen as | ^^(O)) + ^^(O)) ® (- 0 ^( 0 )| + ■ 

This initial state populates the two bands equally with a uniform distribution in k , with 
Ap fc (0) =0. So and fj® described in Eq. (Jb) have no contributions to adiabatic 
pumping. Nevertheless, the effect of IBC on the transport is nonzero, with the agree¬ 
ment between theory and numerics presented in Fig. [2j)a) . It is seen that as dephasing 
strengthens, Q increases first and then decreases, reflecting a competition between dephas¬ 
ing and IBC. I 11 particular, the IBC induced transport is significant even when the dephas¬ 
ing rate 7 is comparable to the characteristic scale of the system’s band gap. Turning to 


a second density matrix as the initial condition, i.e., 

r(“) 


VT6|tf } (0)) + V0Ae ik \^\0)) k 


The term for this case turns out to make no con- 


\/0^(V4 1} (0)| + V0Ae~ ik (^\ 0)| 

tribution. Both theoretical and numerical results for this example are shown in Fig. |2](b). 























o* 
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FIG. 2: (color online) Number of pumped particles Q vs 7 . In (a), the initial state equally populates 
the two bands with IBC and in (b), the initial state unequally populates the two bands with IBC. 
See the main text for details of the initial states used. The system Hamiltonian is described in 
Eq. with 5 = —0.5 in (a) 5 = —1.6 in (b). The adiabatic sweeping speed is chosen to be 

v = 10 -3 in our numerics. In (b), (blue) dotted line, (green) dashed line, and (red) dot-dashed line 
represent the respective contribution from and in Eq. (4); the black line represents 

the total pumping as compared with the numerical results represented by discrete points. 


There, we have separately plotted the contributions associated with , f^ c J, and 
Though the overall dependence of Q on 7 shows one valley only, it is seen that each of the 
three terms responds to dephasing in different manners. That is, it is necessary to know all 
these terms in order to better understand the overall pumping. 

Concluding remarks. Quantum adiabatic pumping with rather general nonequilibirum 
initial states and in the presence of dephasing is contributed by two components of differ¬ 
ent nature, with one of them depending on initial Bloch state populations and the other 
determined by interband coherence. Though our theory uses a very simple pure-dephasing 
Lindblad master equation, the explicit results obtained here should be useful towards under¬ 
standing the dynamics of quantum pumping in more realistic situations. By using different 
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initial states or different pumping protocols, it is possible to isolate each of the two pumping 
components. We highlight the interband coherence effect, as it persists well in the presence 
of dephasing. Note that the interband coherence effect on pumping [fjfj in Eq. (4)] is de¬ 
termined by at the start, so its magnitude can be extensively manipulated by varying 

the switch-on behavior of a pumping protocol. This promises an interesting and relatively 
robust means to coherently control the pumping dynamics albeit dephasing. 
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Appendix 

This Appendix has three sections. In Appendix A, we present detailed derivations of 
Eq. (1) and Eq. (2) in the main text. This is followed by Appendix B, where numerical 
results are described and compared with theoretical results. The most important section is 
perhaps Appendix C, which gives a detailed derivation of the four subterms of fk s used in 
the main text. Whenever possible, we use the same notation as in the main text. 


Appendix A: Details about p m j(s) and Pjj(s) 


We first present again the Lindblad master equation used in the main text, 


^pO) = -i [H(s)p(s) — p(s)H(s)} 


7 \ A(s)p(s)A(s) - - [A 2 (s)p(s) +p(s)A 2 (s)] \ . 


(Al) 


Projecting both sides of this master equation onto |m(s)) and (j(s)|, the instantaneous 
eigenstates of H(s), one obtains 

f <i(s)lj^p(s)|m(s)> = -i(j{s)\H(s)p(s)-p(s)H(s)\m(s)) 

+ 7 - \ P 2 (s)pM + pMVw] |m(»)) (A2) 


Note that if j = m, then all the terms on the right hand side (RHS) of Eq. (A2) will 


disappear. As introduced in the main text, p(s) is now expressed in the representation of 
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instantaneous energy eigenstates of H(s), with p(s) = Y^i n P ln ( s )\K s ))( n ( s )\- Plugging this 
expansion of p(s) into Eq. (A2), we find the diagonal density matrix elements 


Pjj( s ) 

m^j L 


Pmj(s)(j(s)| — \m{s)) +c.c. 


(A3) 


Here means the derivative with respect to the system parameter s. For an adiabatic 
process starting at s = s 0 and ending at s — .sq, the final population on level j is then given 
by 


/ Si 

, 

■rn^j 3 t 


Pmj(s)(j(s)| — |m(s)) + C.C. 


(A4) 


In a similar manner, letting m ^ j in Eq. (A2), we hnd the off-diagonal density matrix 
elements, 


Pmj (s) V 


{m(s)\p(s)\j(s )) 


i [E m (s ) - Ej(s)} - l [A m (s) - A,-(s)] 


2 • 


The numerator in Eq. (A5) can be further rewritten as 


(AS) 


(m(s)|p(s)|j(s)) = Pmj(s) + [ pjj(s ) - Pmm(s)} (m(s)\ — \j(s)) 


E 

n ¥ z m,j 


Pnj(s)(m(s)|4|„(s)) - Pmn (s)(„( s )|il|j( s )) 


(A6) 


Reexpressing the off-diagonal density matrix elements p n j(s) and p mn (s) in the second line 


of Eq. (A6) using Eq. (A5), and then inserting Eq. (k\6b back into Eq. (AS), we hnd that, 


to the first order of v, the terms in the second line of Eq. (A6) can be dropped. So up to 
the first order of v, we finally have 


/ \ _ ^ Pm.j ( s ) + Pmm( s )} (m(s) | ds \j (s)) 

PmAS V -i [EM - E 3 {s )] - 2 [A m (s) - A&)] 


2 • 


(AT) 


Further plugging Eq. (A7) into Eq. (A4), we arrive at 


P„M-Pd») = E f { A-f).°' M , l ^ 1 . |m( ; ) > 1 + c.c.>d 8 

" ^Js o \ gmj&hmA^+WmM 

y r s ' j [Pii( S ) -Pmm(s)] 0'(«)l^r i l m ( s )) 

“g Js o I 9mj (s) [7 mj (s) + Wmj (s)] 


(AS) 


+ c.c. > ds , 


where 


'ymj(s) — 2 [A-m('S) Aj(s)] , 9mj[s) — -^ m (s) Eji^s). 


(A9) 
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To proceed further, we perform an integration by parts over p m j(s), the first term on the 


RHS of Eq. (A 8 ). That term then becomes 


E 

m^j 





+ c.c. 


['ymj(s') “I - ^9mji^S)\ 


s=s i 


S=Sq 


(A10) 


+ c.c. > ds. 


Note that for a nonzero dephasing rate 7 , all the off-diagonal elements p m j(s) would have 
decayed to zero at the end of a slow pumping cycle. Thus, we can ignore the term propor¬ 


tional to p m j(s 1 ) in Eq. (A10). Using again Eq. (A7), one sees that the term on the second 
line in Eq. (A10) is at least of order 0(v 2 ). Thus this term can also be dropped if we only 


consider terms up to the first order of v. With these clarifications, the only term left in 


Eq. (A10) is 


E 

m^j 


9mj(s') T (s)] 


+ c.c. 


(All) 


S=S 0 


The second term in Eq. (A 8 ) can be simplified by replacing Pjj(s) and p mm (s) with their 


zeroth order expressions, namely, Pjj(so) and p mm (so). Together with Eq. (All), we finally 
have 


Pjji s 1 ) Pjj ( s o ) 


m^j 


, dm.j(g)(i(g)| iL S £l |m(s)) 
9mj{s) [Tmj(®) T ^9mj{s)\ 


+ C.C. 


s=s 0 


rsi 7 mj(s) (j(s)\^-\m{s)) 

lto(So) “ ' >rara(So)1 L ”WwhUTUUh 


(A12) 


-ds, 


The expression we show in the main text [Eq. (2) therein] is obtained under the condition 
s(t) = vt. In that case v is a constant and can be pulled out from the above integral in 


Eq. ( |Al2| 


Appendix B: First-order nonadiabatic corrections in two- and three-level systems 


This section is devoted to verifying Eq. (2) in the main text or Eq. (A12) here. We first 
consider a Landau-Zener system (LZS) subject to Lindblad pure dephasing, as described by 


Eq. (A1). The system Hamiltonian is given by H(s) = |( ga(Jx+so z ), where the constant go > 


0 equals the minimum energy gap of the instantaneous Hamiltonian H(s). For convenience 
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we choose A(s) = H(s). As an example, initially (s = so) the two levels of H(sq) are 
assumed to be coherently populated as: 


P(s o) = 


^|^ (1) (s 0 )) + ^ (2) («o)> 


^(V> (I| WI + i» (2) WI 


(Bl) 


where and |-0^(s o )) denote the ground and excited states of H(s) at s = s 0 - In 

our explicit calculations, we choose So = — 1, Si = 1, go = 1. Plugging this initial state into 


Eq. (A12), the nonadiabatic transition probability up to the Erst order of v at the final time 


s = Si then reads: 

A p = p 2 2 {s = 1 ) - P 22 (s = - 1 ) = -v 


a /37 7 f n / 4 cos 4 (6) 

8 (y 2 + 2 ) + 1 2 J_ n/ . 4 y 2 + 4cos 2 (6») 


d0. 


(B2) 


On the RHS of Eq. (B2), the first term is due to initial state coherence and the second 


term is an integral (which can be worked out analytically) over the whole adiabatic process 
weighted by the population difference of the two levels at so = — 1. To numerically verify 


Eq. (B2), we directly evolve the initial state given in Eq. (Bl) using the Lindblad master 


equation [Eq. (Al)] from s 0 = — 1 to Si = 1. The comparison between theory and numerical 
results is shown in Fig. (J3]) , where an excellent agreement is obtained for a wide range of 
the dephasing rates. 

To further check how initial state coherence induces corrections to nonadiabatic popula¬ 
tion transfer, we consider a second initial state: 


POo) = 


A|^< 1 >( S „)) + A|v> ( 2, (V)> 


A{v. (1| (i>„)i + /V 2 >M 


(B3) 


As seen in Eq. (A 12 ) and also stressed in the main text, the switch-on behavior of an 


adiabatic process will be important in correcting nonadiabatic transition probabilities in 
the presence of initial state coherence (hence important to adiabatic pumping). To check 
this we investigate three different adiabatic protocols: (I). s = ut , (11). s = cos (ut), and 
(III), s = 1 — u 2 t 2 . Assuming that the start times are t 0 = —ir/u for case II and t 0 = —y/2/u 
for case III, the sweeping rate of s at the start would be v = u } v = 0 and v = 2a/2 u, 
thus yielding different nonadiabatic corrections. In addition, for this initial state, the two 
levels of H(sq) are equally populated in the beginning, so only the first term on the RHS 
of Eq. ( A12[ ) will give rise to nonadiabatic transitions up to the first order in v. The 
theoretical nonadiabatic transition probabilities for these three protocols are respectively 
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FIG. 3: (color online) Nonadiabatic transition probabilities in LZS under pure dephasing, plotted 
as a function of the dephasing rate 7 . (magenta) dashed line: transition probabilities due to initial 


state coherence [the first term on the RHS of Eq. (B2)], (blue) dotted line: transition probabilities 


due to initial population difference [the second term on the RHS of Eq. (B2)], (green) solid line: 
total transition probabilities in theory, (red) circles: total transition probabilities obtained by 


numerically evolving Eq. (Al). The adiabatic sweeping speed is chosen to be a constant, with 
v = 1(T 3 . 


given by A p = —u ^2 + ^ ? A p = 0, and A p = — 2\/2u 4 ^ +2 ) ■ I 11 Fig. jdj we compare our 
theory with numerical results, with excellent agreement for all the three adiabatic protocols. 

We have also considered a three-level system described by the Hamiltonian H(s) = goS x + 
sS z , where S x and S z are spin-1 operators. In this case, for a constant sweeping rate v of s, 
the nonadiabatic transition probability on the second excited state is given by 


Ap n - —2v [pn(so) - P 22 (so)] D , - vD 2 , 


(B4) 


with 


Di = - A igl 


3 fun 1 


(2-7Wo) cot (0o) 


7 d tan 


3o 


+ 


\A 2 9o+ 4 


\/ 1 2 9o + 4 


+ 


<7o+So 


(B5) 
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FIG. 4: (color online) Nonadiabatic transition probabilities with respect to the dephasing rate 7 
for three different adiabatic protocols. From top to bottom, the adiabatoc protocols are given by 
s = cos (ut), s = ut, and s = 1 — u 2 t 2 , with u = 10 -3 , the starting and ending values of s are 
given by so = —1 and si = 1. Dotted, dashed, solid lines denote theoretical results, while symbols 
denote numerical results. 


and 


Do = 


4v / 2fi'o {Re [pi 2 (s 0 )] 7 9i - 21m [pi 2 (so)]} (^1 + 2s 0 - 1) 


(B 6 ) 


9i + (4 + 7 2 5 , i) V9i + figi + 1) (2s 0 — 1) 
where pi 2 (so) is the initial off-diagonal density matrix element describing the coherence 
between the two excited states and 57 = \J + (1 — 2 s 0 ) 2 . 

Figure [5] shows again the agreement between theory and numerics for this case, for a 
wide range of the dephasing rate 7 and for various initial states. Comparing the two cases 
shown in the main panel of Fig. [5j we indeed see the impact of initial state coherence on 
nonadiabatic transition probabilities, albeit dephasing. 
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FIG. 5: (color online) Nonadiabatic transition probabilities A pu vs. the dephasing rate 7 . 
For two types of initial states, namely, a superposition state of three instantaneous eigenstates 
of H(s 0 ) (hence with initial state coherence) and a mixed state without initial state coher¬ 
ence. In the former case, the initial wavefunction of the system is assumed to be |^(so)) = 
\/o( 8 |l(so)) + V / OA| 2 (so)) + \/o7T|3(so))- In the latter case, the initial density matrix is given by 
p(«o) = 0.8|l(so))(l(so)| + 0.1|2(so))(2(so)| + 0.1|3(so))(3(so)|. In the calculations we set go = 1 , 
v = 10” 3 , so = — 1 and si = 1. Inset shows the nonadiabatic transition probabilities if the initial 
state is prepared in a pure state |i/>) = 11), namely, the second excited state. In all the cases our 
theory agrees with numerics. 

Appendix C: Detailed Derivations of f ks 

As defined in the main text, the number of pumped particles over one adiabatic cycle is 
given by Q = A dk ds f ks , with f ks defined by 

fks = Tr [n _1 p(s)u fe (s)] . (Cl) 
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Here v k (s) = dH j k s ^ is the group velocity operator. Without loss of generality we assume the 
adiabatic pumping protocol starts from s 0 = 0 and ends with si = 27r. Using 


fks = v 1 ^2p jj (s)(j(s)\v k (s)\j(s))+v 1 J2p mj (s)(j(s)\v k (s)\m(s)) (C2) 


], m 
m^j 


as well as the expressions for Pjj(s) and p m j(s ) in Eqs. (A7) and (A12), one obtains 


Ik, = v 1 ~ 2 E M°) ~ P' 




m^j 


dEj(s) r |(j(g / )| li ^|m(5 / ))| 2 7 m j(5 / ) 
dk Jo g 2 mj {s') [l 2 mj {s') + g 2 mj (s')\ 




m^j 


9mj(s ) \pfmj(s ) A )] 


s'=0 


_ | m{c)) A Pmj (g) + [( m ( g )l£l n ( g ))dnj(g) - Pmn(g)(w(s)|^|j(s))] ^ 


m^j 


dk 


'ymj(s') A PJni j (&) 


As emphasized in the main text, all the quantities in the above equation should be under¬ 
stood as functions of quasimomentum k, though this dependence is not explicitly spelled 


out. Let us first focus on the first term on the RHS of Eq. (C3). In principle this term (in¬ 


versely proprotional to v) can have a large contribution to the pumped number of particles. 
To highlight other contributions due to nonadiabatic transitions, we assume that both the 
initial state and the spectrum of H k {s) and A k (s) are even functions of k. Under this as¬ 
sumption, the first term upon integration over k will have no contribution to the pumping. 


For the same reason, the second term on the RHS of Eq. (C3) will also vanish. Under this 
simplification, f ks reduces to 


f _ f( c ) I Aabd) 
J ks J ks ' J ks ’ 


(C4) 


with 


and 


r(abd) _ 

J ks 


r( c ) __ 
J ks 


m^j 


SEM 

dk 


A Pmj ) 


9mj(s )[T / m i j('S ) A )] 


(C5) 


s '=0 


-E 

m^j 


(j ( s ) | dg ( g ) | m ( „)) + En [( m ( g ) I £ K S )>Pnj(s) - p mn (s) (n(s) \ £ | j(s))] 


dk 


'Hmj(s') A ifjmj (s) 


(C6) 
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Next we focus on the first term on the RHS of Eq. (06), which contains j-p mj (s). Because 
fks will be under an integration upon s to give the number of pumped particles Q, we can 
consider an integration over s by parts here. Then this term (under the integration over s ) 
will result in two terms, 

ywi^hw) 


-E 

j,m ' 
m^j 


O'(s)l^jpMs)). 




Pmj (s)ds — ^ ^ 


m^j 


Pmj (‘ 


s=2n 


s =0 


/»27T 

'y 1 / Pmj 


iQmj (-S) 


m/j 


(C7) 


ds. 


ds 7 mj(s) + ig mj (s) 

Since a nonzero dephasing rate 7 leads to an exponential decay of all off-diagonal density 


matrix elements p mj (s), the term proportional to p m j{s ) in Eq. (C7) is negligibly small at 
the end of a slow pumping cycle (s = 2n). Thus the only contribution of the first term on 


the RHS of Eq. (C7) to the pumping current is due to: 

(j( s ) 7r 1 l m ( s )> 


E 

j,m 

m^j 


Pmj ( s ) 


r )mj{^S) H - ^Qmj (<s) 


(C 8 ) 


s=0 


The second term in Eq. (C7) still involve an integral over s. The integration of the off- 


diagonal density matrix element p m j(s ) (oscillating with s while decaying) over s indicates 


that it is at least of the order of v, hence negligible as compared to the first term in Eq. (C7). 


Alternatively, one may plug our earlier expression for p m j(s ) in Eq. (A7), only to find that 


the second term in Eq. (C7) becomes 


7 mj(s) + igmj(s) 


d (jPOI^prdraQ)) 

ds 7 mj{s)+ig m j(s) 


ds, (C9) 


which is indeed of the order of v. As such, depicted in Eq. (C 6 ) equivalently (upon 


integration of s over one pumping cycle) contains the following contributions that have 
u-independent terms, 


n(abd) 
'ks 


where 


Ad) _ 

Jks 2 tt^ 


Pmj (* 


J,m l 
m^j 


Ad) Aab) 

J ks ' J ks "> 

(CIO) 

j0'(«) 1 dfcV(s)) 

7 mj(s) T 'mj(s) s= o 

(Cll) 
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and 


fl? = -£owi 


3,m 

m^j 


d H(s) 
dk 


|m(s)) 


En [( m ( S )l£l n ( g ))^( g ) -Pmn(3)(w(g)|^|j(g))] 
T / mj('S) “1“ djnij ( S ) 


• (C12) 


As we have explained earlier, the contributions of those terms involving the off-diagonal 
density matrix elements to the pumping are at least of order O(v) after integrating over s. 


Ignoring these higher-order effects, Eq. (C12) reduces to 

kijM - 


/r=-E«t>i 


3,m 

m^j 


d H{s) 
dk 


m(s))- 


d~ d] r nj (S ) 


(C13) 


With the help of g dk ^" (A ^ = (j(s)|^|m(s)), this above expression for becomes 


fl? = -E 


3,m 

m^j 


\pjj( s ) ~ Pmm(s)} (m(s)|^|j(s))(j(g)|^|?n(g)) 

r m j(s) + i 


(C14) 


where T m j(s) = ^j is dehned in the main text. To proceed further, we shift the derivative 
of k to act on state | j(s)) and exchange m and j in the summation. With these manipula¬ 


tions, Eq. (C14) is converted to 

fks ] = E^( s ) 


3 ,rn 

m^j 


(^'(s)|m(s))(m(s)|£j(s)) + (£j(s)\m(s))(m(s)\±](s)) 


f ^mj (s) T i 




(CIS) 


Recombining the two terms in Eq. (C15) and using the zeroth order expression for Pjj(s) (be¬ 
cause this only introduces an error proportional to v to the pumping, which can be neglected 
in a slow pumping protocol), we arrive at 


f ( ab ) __ fW I f(°) 

J ks J ks ' Jks J 


0 ) i f(P) 


(C16) 


where 


and 


r (a) _ 0 v^ /m W s ) Re [(^i(s)|m(s))(m(s)|^i(s))] 

hs r2-(s) + l 

7 ,m " L J 


Ab) 

J ks 


0 ) 

m^j 


Im [(abX g )M g ))( m ( g )l£7'(g))] 

W + i 


(C17) 


(CIS) 


Combining all these results together, we finally have fks = + fks + fjfj + flt\ which 

is precisely the result used in the main text. It should be emphasized that all the four 
v-independent subterms of fk s are derived based on nonadiabatic corrections captured to 
the first order of v. 
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